Dynamics of nuclear single-particle structure in covariant theory of particle-vibration 

coupling: from light to superheavy nuclei 
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The impact of particle-vibration coupling and polarization effects due to deformation and time- 
odd mean fields on single-particle spectra is studied systematically in doubly magic nuclei from 
low mass ^''Ni up to superheavy ones. Particle-vibration coupling is treated fully self-consistently 
within the framework of relativistic particle-vibration coupling model. Polarization effects due to 
deformation and time-odd mean field induced by odd particle are computed within covariant den- 
sity functional theory. It has been found that among these contributions the coupling to vibrations 
makes a major impact on the single-particle structure. The impact of particle-vibration coupling 
and polarization effects on calculated single-particle spectra, the size of the shell gaps, the spin-orbit 
splittings and the energy splittings in pseudospin doublets is discussed in detail; these physical ob- 
servables are compared with experiment. Particle- vibration coupling has to be taken into account 
when model calculations are compared with experiment since this coupling is responsible for ob- 
served fragmentation of experimental levels; experimental spectroscopic factors are reasonably well 
described in model calculations. 

PACS numbers: 21.10.Pc,21.10.Ro,21.60.Jz,27.90.+b 
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I. INTRODUCTION 

The covariant density functional theory (CDFT) [3-11] 
is one of standard tools of nuclear theory which offers 
considerable potential for further development. Built on 
Lorentz covariance and the Dirac equation, it provides 
a natural incorporation of spin degrees of freedom [l|, 3 
and an accurate description of spin-orbit splittings 2| 
(see also Fig. 2 in Ref. Q), which has an essential in- 
fluence on the underlying shell structure. Note that the 
spin-orbit interaction is a relativistic effect, which arises 
naturally in the CDFT theory. Lorentz covariance of the 
CDFT equations leads to the fact that time-odd mean 
fields of this theory are determined as spatial components 
of Lorentz vectors and therefore coupled with the same 
constants as time-like components [5[ which are fitted 
to ground state properties of finite nuclei. In addition, 
pseudo-spin symmetry finds a natural explanation in the 
relativistic framework Q. CDFT in its different incar- 
nations both on the mean field and beyond mean field 
levels provides succesful description of many properties 
of ground state and excited configurations in nuclei Q. 

However, the majority of applications of CDFT have 
been focused on collective properties of nuclei. There are 
only few features on nuclear systems, strongly dependent 
on single-particle degrees of freedom, which have been 
addressed in the CDFT studies on the mean field level 
and compared with experiment. These are the single- 
particle properties in spherical and deformed nuclei, spin- 
orbit splittings in spherical nuclei, magnetic moments in 
the A±l neighbours of the ^^O nucleus {A = 16) [?[, and 
alignment properties of single-particle orbitals in rotating 



nuclei [8j\ 

The calculated single-particle spectra of spherical nu- 
clei are frequently compared with experiment ( see, for 
example, Refs. [l|, 3). However, only in Refs. 0, [H de- 
formation polarization effects and time-odd mean fields 
are included in the calculations which makes them more 
realistic. The spin-orbit splittings, as extracted from 
the single-particle energies obtained in spherical CDFT 
calculations, reproduce the experimental spin-orbit split- 
tings fairly well, although there are deviations up to 20% 
of absolute value of the splitting (see Figs. 11-12 and 
Table 4 in Ref. [H and Fig. 2 and Sect. IVB in Ref. 0) 
which only weakly depend on the RMF parametrization. 
Note that the comparison of Ref. [3l includes the un- 
certainty of the spin-orbit splittings due to deformation 
polarization effects and time-odd mean fields as they are 
found in Ref. [l^. Ref. Q clearly shows that spin-orbit 
splittings are better described in CDFT as compared 
with Skyrme energy density functional theory (SEDFT) 
despite the fact that no single-particle information (con- 
trary to SEDFT) has been used in the fit of the RMF 



^ Effective alignments of compared rotational bands, which sensi- 
tively depend on both the alignment properties of single-particle 
orbital by which two bands differ and polarization effects induced 
by the particle in this orbital, are in average better reproduced 
in the CDFT calculations than in the cranked Nilsson-Strutinsky 
calculations based on phenomenological Nilsson potential, see 
comparisons presented in Refs. [9l 4lll . This is despite the fact 
that no single-particle information has been used in the fit of the 
CDFT parameters, while the parameters of the Nilsson potential 
are fitted to experimental single-particle energies (see Ref. Il2ll ). 
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FIG. 1: (Color online) Single-particle strength distributions of the neutron lg9/2 state in "^sn calculated within the PVC 
model with different numbers of the phonon modes. Red bars show the initial mean field lg9/2 neutron state. 



Lagrangian. 

These comparisons, however, do not reveal the accu- 
racy of the description of the single-particle states be- 
cause the particle- vibration coupling (PVC), which can 
affect considerably the energies of single-particle states in 
odd-mass nuclei has been neglected. The mod- 

ification of the quasiparticle states by par ticle- vibration 
coupling is weaker in deformed nuclei [l9l [20 l| since the 
surface vibrations are more fragmented (less collective) 
than in spherical nuclei [l^, [2l|. As a consequence, the 
corrections to the energies of quasiparticle states in odd 
nuclei due to particle- vibration coupling are expected to 
be less state-dependent in deformed nuclei. Hence the 
comparison between experimental and calculated ener- 
gies of single-particle states is expected to be less am- 
biguous in deformed nuclei as compared with spherical 
ones [H, at least at low excitation energies, where 
vibrational admixtures in the wave functions are small. 
The analysis of the energies of one-quasiparticle deformed 
states in actinide region, performed within the relativis- 
tic Hartree-Bogoliubov (RHB) approach in Refs. [2^[23|. 
reveals that while the majority of the states are described 
with an accuracy better than 0.5 MeV, there are a num- 
ber of states which deviate from experiment by as much 
as 1 MeV. 



The question arises of which features of the DFT are 
most decisive for the single-particle structure. The most 
crucial ingredients in this respect are 

• the effective nucleon mass and its radial dependence 
which determines the level density near the Fermi 
surface. It is well known that the one- (quasi-) par- 
ticle spectra calculated on the mean field level are 
less dense than in experiment both at spherical [l^ 
[iTj and deformed ^2,23] shapes. The average level 
density of the single-particle states on the mean 
field level is related to the effective mass (Lorentz 
mass in the notation of Ref. [IJ] for the case of 
CDFT theory) of the nucleons at the Fermi surface 
m\{kF)/rn. The CDFT is characterized by a low 
effective mass m*^{kF)/m « 0.66 [1]. As a conse- 
quence, this low effective mass of CDFT has a big- 
ger impact on the single-particle spectra than the 
effective masses of non-relativistic DFT which are 
typically larger. Note that at spherical shape the 
inclusion of particle- vibrational coupling brings the 
calculated level density closer to experimental one 
characterized by an effective mass m*{kp)/m w 1.0 

• the spin- orbit potential which determines the ener- 
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FIG. 2: (Color online) Spectra of ^^Ni and its neighboring odd nuclei. Column 'sph' shows the single-particle spectra obtained 
in spherical RMF calculations of ^^Ni. Column 'sph+PVC shows the spectra obtained in spherical calculations within the 
PVC model. Columns "def", "def+TO", "hybrid" and "exp" show one-nucleon separation energies defined according to Eqs. 
([T]) and ((2]). Column "def" is based on the results of triaxial CRMF calculations with no TO mean fields. These fields are 
included in the calculations the results of which are shown in column "def+TO". The corrections due to PVC are added in 
column "hybrid" . The experimental single-particle energies are displayed in column "exp" ; they are based on the data of Refs. 
[5^ (masses of ground states), [sl] (^^Ni), [s!] ('^^Co), [30,|11| (^^Ni) and fsO] ('"^Cu). In order to distinguish overlapping levels, 
orange and then maroon colors are used for the levels in addition to their standard color used in a given column. See text for 
more details. 



getic distance of the spin-orbit partners. It is well 
known that CDFT describes rather well spin-orbit 
splittings [4]. However, these results have been ob- 
tained on the mean field level, and it is necessary 
to test how they will be affected by the inclusion of 
particle-vibrational coupling. In CDFT theory the 
Dirac effective mass m*jj{kF)/rn is closely related 
to the effective spin-orbit single-nucleon potential, 
and empirical energy spacings between spin-orbit 
partner states in finite nuclei determine a rela- 
tively narrow interval of allowed values: 0.57 < 
m'^{kF)/m < 0.61 [2^. This requirement restricts 
considerably the possible parametrizations of the 
RMF Lagrangian on the mean field level, and thus 
it would be interesting to see whether the inclusion 
of PVC will loosen this requirement and thus allow 
wider range of effective masses in the CDFT. 



• the density dependence of the effective potential and 
effective mass which has an influence on the relative 
position of the states with different orbital angular 
momentum I. The investigation of Ref. [13], cov- 
ering spherical odd nuclei bordering doubly magic 
nuclei clearly showed that the relative placement of 
the states with different orbital angular momenta 
/ is not so well reproduced on the mean field level 
in CDFT. The problems occur predominantly in 
connection with the states of high I which hints 
that the surface profile of the mean-field and ki- 
netic terms are involved. Microscopic consideration 
in many-body theory indicates that the effective 
mass has a pronounced surface profile which is in- 
sufficiently parametrized in the present mean-field 
models flG]. Thus, it is important to see whether 
the inclusion of PVC will improve the description 
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of the relative positions of the states with different 
orbital angular momentum I. 

In Ref. a relativistic particle- vibration coupling 

model has been developed and the calculations based on 
the NL3 [2^ parametrization of the CDFT have been 
performed for the doubly magic ^°*Pb nucleus. In this 
work we present a more systematic investigation within 
this model including the calculations covering the nu- 
clei from light ^^Ni and medium-mass ^"^Sn, ^■^^Sn to 
heavy ^"^Pb and superheavy ^^^120 doubly magic nu- 
clei using improved (as compared with NL3) NL3* [27| 
parametrization of the CDFT. This investigation covers 
neutron-deficient near proton drip-line nuclei (^^Ni and 
^"^^Sn), neutron-rich ^^^Sn nucleus, and the ^"^Pb nu- 
cleus located at the beta-stability line. In addition, the 
polarization effects in odd-mass nuclei due to deforma- 
tion and time-odd mean fields are treated by means of 
cranked relativistic mean field (CRMF) calculations. 

The article is organized as follows. The method of com- 
parison of experimental and theoretical single-particle 
levels is discussed in Sec. |lll Section [ITT] presents rela- 
tivistic particle- vibration model and cranked relativistic 
mean field theory and their details related to the study of 
single-particle states. The comparison of model calcula- 
tions with experiment is presented in Sec. IIVI The accu- 
racy of the description of single-particle spectra and spec- 
troscopic factors, the impact of particle- vibration cou- 
pling on shell gaps, spin-orbit splittings and pseudo-spin 
doublets as well as the role of polarization effects in odd- 
mass nuclei due to deformation and time-odd mean fields 
are discussed in detail. Finally, Sec. |V] reports the main 
conclusions of our work. 



II. EXPERIMENTAL VERSUS THEORETICAL 
SINGLE-PARTICLE ENERGIES 

The experimental information about single-particle 
level structure of even-even doubly magic nucleus with 
the proton number Zq and neutron number A^o originates 
from the odd-mass {Zq ± 1) and (Nq ± 1) nuclear neigh- 
bors. In odd-mass nuclei, the polarization effects due 
to deformation and time-odd (TO) mean fields^ induced 
by odd particle play an important role (see Sec. I IV Bp . 
In addition, the particle- vibration coupling modifies the 
energies and the structure of the levels (Sees. IIVCI and 
IIV Ep : neither of them remains purely single-particle in 
nature. In order to describe all these effects we use a hy- 
brid model (labelled as "hybrid" below) , in which the po- 
larization effects due to deformation and TO mean fields 



^ The name nuclear magnetism is frequently used in the CDFT 
to describe the effect of time-odd mean fields induced by mag- 
netic potential of the Dirac equation since this potential has the 
structure of a magnetic field [3, |2S(I . 



are treated on the mean field level within cranked rel- 
ativistic mean field (CRMF) approach (see Sect. IIII Ap . 
while the corrections due to the PVC are treated within 
relativistic particle-vibration coupling model (see Sect. 
IIII Bp . The need for a such hybrid model is dictated by 
the fact that the PVC model is not variational in nature 
and thus it does not allow to calculate polarization ef- 
fects due to odd particle. Since the polarization effects 
and the corrections due to the PVC have a different ori- 
gin, they are treated as additive quantities in the hybrid 
model which is a reasonable approximation. The details 
of the hybrid model are discussed in Sec. IIIII 

The experimental and calculated energies of the parti- 
cle [e (particle)] and hole [£(hole)] states closest to the 
Fermi level are determined from the difference of the 
binding energies of the core [i?(core)] and the corre- 
sponding adjacent odd nuclei [i3(core + nucleon) and 
i3(core — nucleon)] as [isl [29l| 

e (particle) — B{core) — B{core + nucleon) (1) 

and 

e(hole) = i3(core — nucleon) — i3(core). (2) 

These quantities correspond to one-particle removal en- 
ergies. We will call them as single-particle energies for 
simplicity of discussion. Doubly magic spherical nuclei 
^^Ni, ^°°'-^^^Sn, and ^°*Pb are used as the cores in our 
analysis. Their ground states are not affected by PVC. 
The experimental energies of particles and holes defined 
according to these equations include the polarization ef- 
fects due to deformation and TO mean fields induced by 
odd particle or hole as well as the energy corrections due 
to particle-vibration coupling in odd-mass nuclei. Thus, 
their comparison with the results of "hybrid" calculations 
is straightforward. 

On the other hand, there are some unresolved ques- 
tions when these energies are compared with pure mean 
field calculations as it is done frequently in the litera- 
ture. The problem is related to the fact that experi- 
mental single-particle levels and their energies and struc- 
ture are affected by the particle- vibration coupling which 
is neglected on the mean field level. Each mean field 
state k with the energy Sk is split into many levels 
due to particle-vibration coupling, so the single-particle 
strength is fragmented over many levels. In the diagonal 
approximation for the nucleonic self-energy, these levels 
have the same quantum numbers as the original mean 
field state k, but different energies and spectroscopic 
factors S^. Unlike the Hartree or Hartree-Fock approxi- 
mations without pairing where the single-particle states 
are either fully occupied or empty, in the PVC model 
their occupation probabilities are the real numbers be- 
tween zero and one, so that sum rule 

E^^— 1 (3) 

is satisfied. For the states in the vicinity of the Fermi 
surface one usually obtains one dominant level with 
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0.5 < 5^ < 1.0 and many other levels with small S^. 
For the states far from the Fermi surface one observes 
very strong splitting over many levels with much smaller 
and comparable spectroscopic factors, see the detailed 
analysis in Ref. [13] for ^°^Pb. 

This discussion clearly indicates that some procedure 
has to be employed in order to extract pure (or "bare") 
single-particle energies from experimental data which 
have to be compared with mean field single-particle en- 
ergies £k- Two such procedures have been discussed in 
the literature. 

In the first procedure, the energy of the dominant frag- 
ment is corrected by the energy shift due to particle- 
vibration coupling leading to a "bare" single-particle en- 
ergy. Such energies of the doubly-magic nuclei are ex- 
tracted in some publications, see, for example, Ref. [30| . 
However, this procedure relies on the choice of particle- 
vibration coupling model. However, different PVC mod- 
els lead to different energy shifts (see Table 1 in Ref. 



3 If ) . In addition, such procedure frequently neglects po- 
larization effects due to deformation and TO mean fields 
which are again model dependent. For example, there is 
no time-odd mean fields in the phenomenological models 
based on the Nilsson or Woods-Saxon potentials. 

In the second procedure, the center of gravity of the 
fragmented levels with a given quantum number j'^ [j 
is the total angular momentum and tt is a parity) and 
energies is obtained via weighted-average [32, ISj] pro- 
cedure 



(4) 



This energy is then associated with a "bare" single- 
particle energy. This expression is fulfilled exactly in the 
PVC model discussed in Sec. lIIIBl and it is satisfied with 
high accuracy in this model numerical implementation. 
The normalization in this procedure (which is omitted in 
some publications) is used to minimize the uncertainties 
due to large experimental errors on spectroscopic fac- 
tors, since these errors lead to the fact that the sum rule 
'S'fc = 1 is poorly fulfilled in the majority of experi- 
ments. The application of this procedure requires that 
all i/-levels, to which the single-particle state k is frag- 
mented, are identified in the experiment and that spec- 
troscopic factors are measured accurately. These condi- 
tions are definetely not satisfied in the ^^Ni and ^"^'-'^'^^Sn 
nuclei. 

The situation is somewhat better in ^"^Pb, for which 
the e^''"" values have been extracted in Ref. [s^ for a 
number of levels. However, even in this case the accu- 
racy of the definition of s^^""^ is not known because of the 
number of the reasons listed below. The absolute values 
of spectroscopic factors are characterized by large ambi- 
guities and depend strongly on the reaction employed in 
experiment and the reaction model used in the analysis 
[33I [3^ . The spectroscopic factors in odd mass nuclei 
have been obtained in different reactions: for example, 



in the (d,3He) reaction for 207t1 [3|, (a, t) for ^oSBi 
ill, (^He, a) m and {d,t) [H for 207pb. Even in the 
same nucleus different reactions give different spectro- 
scopic factors (see, for example, Table I in Ref. 38] and 
Table Hill in the current manuscript). In addition, the 
spins and parities of many high-lying fragments are un- 
certain (see, for example. Table 1 in Ref. [ssl]) and/or 
only the part of the sum rule strength has been observed 
in experiment (see, for example, Ref. (s^). As a conse- 
quence, it is not clear whether the complete set of the v- 
Icvels, to which the single-particle state k is fragmented, 
and their spectroscopic factors have been identified in 
experiment. 

In the light of this discussion it is clear that "bare" 
single-particle energies cannot be defined with control- 
lable precision in the nuclei under study. As a result, 
we compare observed experimental levels directly with 
mean field calculations. Such an approach is frequently 
used in the literature (see, for example, Refs. [2, 13, 29]). 
It has clear limitations, but at least it allows to see (i) 
whether the inclusion of particle-vibration coupling im- 
proves the description of specific physical observables and 
(ii) whether the conclusions reached earlier on the mean 
field level are valid or not. 



III. THEORETICAL FRAMEWORK 

A. Mean field level: polarization effects 

Mean field results related to polarization effects due 
to deformation and TO mean fields in odd-mass nuclei 
have been obtained usin g cr anked relativistic mean field 
(CRMF) theory [H, [H, EO] . Ahhough this theory has 
initially been developed for the description of rotating 
nuclei, it is also able to describe the nuclear systems with 
broken time-reversal symmetry in intrinsic frame at no 
rotation in which TO mean fields play an important role 
(see Ref. Q). 

The application of the CRMF theory has the following 
advantages. First, the CRMF computer code is formu- 
lated in the signature basis. As a result, the breaking of 
Kramer's degeneracy of the single-particle states in odd- 
mass nuclei is taken into account in a fully self-consistent 
way. This is important for an accurate description of 
TO mean fields (Ref. [5]). These fields have consider- 
able impact on the moments of inertia in rotating nuclei 
(4l| and the fact that the properties of such nuclei are 
well described in the CRMF code (see Refs. 0, HlJ adds 
confidence to proper description of TO mean fields in 
non-rotating systems [i]. Second, the CRMF computer 
code is formulated in three-dimensional cartesian coordi- 
nates which allows to describe not only axially deformed 
but also triaxial nuclear shapes. This is important for a 
proper description of polarization effects due to deforma- 
tion. 

It is well known fact that the description of pairing cor- 
relations in doubly magic nuclei and their close neighbor- 
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FIG. 3: (Color online) The same as in Fig. [2]but for the spectra of ^""Sn. Extrapolated "experimental" single-particle energies 
of proton and neutron spherical subshells are taken from Ref. [s^ . 



hood is notoriously difficult and model dependent (Ref. 
(i^l). Tlie presence of large sliell gaps leads to a pairing 
collapse in the models with the treatment of pairing on 
the BCS or HFB [Hartree-Fock-Bogoliubov] levels. This 
is seen, for example, in non-relativistic HFB calculations 
with the Gogny DIS force of io"^i32gj^ -^^ ^LeL [42] . The 
calculations within the relativistic Hartree-Bogoliubov 
(RHB) framework of Ref. jl^ with the same force for 
pairing also show pairing collapse in the ground states of 
^^Ni, ioo,i32gj^ j^^j 208p|3^ rpj^g pairing correlations are 

restored when the particle number projection is imple- 
mented. However, there is a substantial difference be- 
tween different approaches (such as Lipkin-Nogami (LN) 
method of approximate particle number projection or 
projection after variation (PAV) approach) and variation 
after projection (VAP) approach; the latter represents 
the best approach for the treatment of pairing correla- 
tions [l^. The RHB-fLN calculations with the Gogny 
DIS force for pairing show very weak pairing in the nu- 
clei of interest. The additional binding due to pairing 
(relatively to the results with no pairing) is only 0.333 
MeV, 0.284 MeV, 0.493 MeV and 0.557 MeV in ^^Ni, 
ioo,i32gj^ and ^°^Pb, respectively. The blocking of single- 



particle orbital in odd-mass nucleus leads to the weak- 
ening of pairing correlations and to frequent unphysical 
pairing collapse even in the RHB-I-LN calculations. The 
physical quantities, such as given by Ecj^s. ([1]) and ([5]), 
are relative in nature. Thus, we believe that in the case 
of weak pairing they are sufficiently well described in un- 
paired calculations. As a consequence, we neglect pairing 
correlations in the CRMF calculations. Note that they 
are also neglected in the PVC calculations. 



Similar to Ref. [1^, binding energies of odd- A nu- 
clei were calculated by blocking the {rir -t- 1, j'^, j) state 
(amongst the (n^ -I- 1, j'^, ra) states) with the largest pro- 
jection m = j onto the symmetry axis. Here, the and 
j are the single-nucleon radial and total angular momen- 
tum quantum numbers. The spherical subshell label of 
this state is used to label the configuration of odd mass 
nucleus since the deformations (induced by odd particle 
or hole) of immediate odd-A neighbors of doubly magic 
nuclei are rather small. 
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B. Beyond the covariant density functional theory: 
single-particle spectra of nuclei 

1. Particle-vibration coupling model of the time dependence 
of the nucleonic self-energy 

The mean field approach, being a very useful and con- 
venient theoretical tool to describe finite nuclei, is, how- 
ever, only a static and local approximation to the many- 
body problem. The non-locality (or, equivalently, mo- 
mentum dependence) is implicitly taken into account by 
an effective mass which is considerably smaller than the 
bare nucleon mass. This enables one to reproduce rea- 
sonably well the general properties of the single-particle 
motion in the nuclear potential well. However, the mean 
field approximations based on the widely used Skyrme, 
Gogny as well as covariant density functionals generate 
too stretched single-particle spectra around the Fermi 
surface and often fail to describe the sequence of the 
single-particle levels. Proper inclusion of the time depen- 
dence in the nucleonic potential requires to go beyond the 
mean field description which is static by construction. 

The underlying physics which is responsible for the 



time dependence is related to many-body correlations, 
i.e. to the part of the hamiltonian which cannot be ex- 
pressed in terms of the one-body density. In the medium 
mass and heavy nuclei the leading order contribution of 
such correlations to the nucleonic self-energy comes from 
coupling of the nucleonic motion to the surface and vol- 
ume vibrational modes (phonons). The most important 
collective vibrational modes are generated by the nuclear 
oscillations of the coherent nature. Taking into account 
the coupling of these modes to the single-particle motion, 
the shell model acquires a dynamic content. These sur- 
face and volume oscillations, especially their low-lying 
modes, modify considerably the picture of the single- 
particle motion in nuclei. The main assumption of the 
quasiparticle-phonon coupling model is that two types of 
elementary degrees of freedom - one-quasiparticle (l-qp) 
and collective ones - are coupled in such a way that con- 
figurations of l-qpig)phonon type with low-lying phonons 
strongly compete with simple l-qp configurations located 
close in energy or, in other words, that quasiparticles can 
emit and absorb phonons with rather high probabilities 

mill. 

Taking the quasiparticle-phonon coupling into account 
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FIG. 5: (Color online) The same as in Fig. [2] but for the spectra of ^''^Pb. 
the data of Refs. f53,E3,[6l. 



The experimental single-particle levels are based on 



as an additional time dependent part of the nucleonic 
self-energy leads to the splitting of each mean-field single- 
particle state into many energy levels. The dominant lev- 
els, i.e. the levels with the largest spectroscopic factors, 
shift, as a rule, towards the Fermi level improving consid- 
erably the agreement with data, see the review of various 
applications in Ref. [IB] and references therein. This re- 
sult is obtained, however, only for the states in about 10 
MeV vicinity of the Fermi energy. For the states lying far 
from the Fermi surface one observes a very strong frag- 
mentation making an extraction of the dominant levels 
impossible. 

Similar general picture has been obtained within the 
particle-vibration coupling extension of the covariant 
density functional theory . Unlike the non-relativistic 
versions of the model, in the extended CDFT the particle- 
phonon coupling self-energy allows intermediate nucle- 
onic propagation through the Dirac sea states with neg- 
ative energies, in addition to the particle and hole states 
in the Fermi sea. The contribution of the Dirac-sea terms 
to the self-energy at the Fermi level are, however, found 
negligibly small due to their large energy denominators. 



2. Approximations 

In the present investigation we retain basically calcu- 
lational scheme of Ref. [17,]. This approach implies lin- 
earized version of particle-vibration coupling model: the 
Dyson equation for single-particle Green function con- 
tains phonon-coupling self-energy with mean-field inter- 
mediate nucleonic propagators and phonon energies and 
vertices computed within relativistic random phase ap- 
proximation (RRPA) [i^. For medium- mass and heavy 
nuclei the RRPA gives a very reasonable description of 
the phonon spectra. Although a more precise description 
of the phonons can be achieved within an approach like 
relativistic time blocking approximation (RTBA) |47j, 
particle-vibration coupling model of Ref. 17} by using 
the RRPA phonons takes into account the major con- 
tribution of the vibrational motion to the single-particle 
spectra. 

The application of the model to the lightest doubly 
magic nuclei "^^O, *"Ca, and "^^Ca requires, however, 
a special consideration because their mean-field single- 
particle level densities are too small to provide a sizable 
configuration mixing around the Fermi surface. For -"^^O 
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and Ca, it is not possible, in principle, to reproduce 
within the RRPA the low-lying phonons with positive 
parity (first positive parity particle-hole excitations are 
at too high energies). In such cases the phonon spectra 
for the nucleonic self-energy should be calculated allow- 
ing an extension for coupling to more complex configu- 
rations; that implies a generalized particle- vibration cou- 
pling model. This will be done in a separate work. In 
''^Ca we have obtained a rather reasonable phonon spec- 
trum, but the strong fragmentation of the hole states 
does not allow extraction of dominant single-particle lev- 
els. 

The phonons of unnatural parities are known to play 
only a marginal role in the nucleonic self-energy, there- 
fore, they are not included into the phonon space. More- 
over, the inclusion of pairing vibrations into the phonon 
space has not been studied yet within the relativistic 
framework. In principle, they may bring some additional 
corrections to the results, although our approximation is 
justified by the results of Refs. SI]) where the frag- 
mentation of the collective modes has been reproduced 
very well without inclusion of the pairing vibrations. 



A rather good convergence of the phonon coupling self- 
energy with extension of the phonon space has been ob- 
tained numerically: with very few exceptions, addition of 
phonon modes with energies above 15 MeV does not af- 
fect the results. Partial contributions from the phonons 
with low angular momenta and low energies have been 
also investigated. An example of such a study is shown 
in Fig. [l]for the neutron lg9/2 state in ^'^^Sn. This typi- 
cal deeply bound hole state is strongly fragmented when 
the PVC is taken into account, see the panel (d) with the 
results obtained with the full phonon set below 15 MeV. 

Panels (a) and (b) show the results obtained when only 
the phonons below 5 and 10 MeV, respectively, are taken 
into account. One can see that in the former case the 
picture is very different from the one displayed in the 
panel (d). In the latter case the results are already close 
to saturation, but there are still two states between —20 
and —18 MeV with the spectroscopic factors larger than 
0.1. Panel (c) shows how the angular momentum cut-off 
weakens the fragmentation effect. Similar conclusion can 
be drawn for other states either far from or around the 
Fermi surface: cut-off on the angular momenta and on 
the energies of phonons has to be made with care. 



C. The numerical details of calculations 

The CRMF equations are solved in the basis of 
an anisotropic three-dimensional harmonic oscillator in 
Cartesian coordinates characterized by the deformation 
parameters (3q — 0.0 and 7 = 0° as well as the oscillator 
frequency huo — 41 A^^^^ McV. The truncation of the 
basis is performed in such a way that all states belonging 
to the shells up to fermionic A^j;-=20 and bosonic Nb—20 
are taken into account in the calculations. Numerical 
analysis indicates that this truncation scheme provides 
sufficient numerical accuracy for the physical quantities 
of interest. Note that the same truncation of basis is 
used in the RRPA and PVC calculations, but they are 
performed at spherical shape. 

The NL3* jl^l parametrization of the RMF Lagrangian 
is used in the CRMF, RRPA and PVC calculations. This 
recently fitted parametrization has been successfully ap- 
plied to the description of binding energies [27j , ground 
state properties of deformed nuclei [4^, fission barri- 
ers do], rotatin g nu clei [27| . giant resonances [13], and 
breathing mode 151]. 



D. Cut-off problem of the phonon basis 

Phonons of the multipolarities 2+, 3^, 4+, 5^, 6+ with 
energies below 15 MeV are included in the model space 
of the PVC calculations. We have neglected only the 
phonon modes with very small transition probabilities, 
less than 5% of the maximal ones for each J'^. The 
phonon energies and their coupling vertices have been 
computed within the self-consistent RRPA. 



IV. RESULTS AND DISCUSSION 
A. Theoretical and experimental spectra 

An example of the comparison of the results of calcu- 
lations and experiment is shown in Fig. [21 Column "sph" 
shows single-particle spectra obtained in spherical RMF 
calculations of even-even ^^Ni nucleus. Column "def" 
displays one-nucleon separation energies [defined accord- 
ing to Eqs. (dJ and ^] obtained in triaxial CRMF cal- 
culations in which only deformation polarization effects 
are taken into account. It is seen that these energies are 
close to the single-particle energies obtained in spheri- 
cal calculations. This shows that deformation polariza- 
tion effects induced by extra particle or hole are rather 
modest. In general, they lead to a slight compression of 
the calculated spectra as compared to those obtained in 
spherical calculations. Note that the magnitude and im- 
pact of deformation polarization effects is similar in our 
calculations and in the calculations of Ref . [l^ . 

The inclusion of time-odd mean fields induces addi- 
tional binding in odd-mass nuclei (see Ref. Q for more 
detail) which is rather modest being around 100-200 keV 
in the nuclei around ^^Ni. This again introduces slight 
compression of the calculated spectra (Fig. (2). One 
should note that additional binding due to TO mean 
fields in odd-mass nuclei (and a subsequent compression 
of the spectra) is substantially larger in the calculations 
of Ref. [I3j. For example, it is around 0.5 MeV in the 
nuclei around ^^'^Pb and reaches 1 — 1.5 MeV in ^^Ni. 
Note that for the latter nucleus additional binding due 
to TO mean fields is only 0.1 - 0.2 MeV in our CRMF 
calculations (see Fig. (2). The reason for this difference is 
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state can be quantified by the energy difference between 
the energies of this state in the columns "sph" and 
"sph+PVC" of Fig. [5] These energy differences treated 
as the corrections due to PVC are then added to one- 
neutron separation energies obtained in triaxial CRMF 
calculations; this leads to column "hybrid". The results 
in this column are compared with experimental data pre- 
sented in column "exp" . 

The single-particle spectra of doubly magic ^°°Sn, 
^■^^Sn and ^"^^Pb nuclei, obtained in the same way as Fig. 
[21 are displayed in Figs. [3l |4] and [5) They show simi- 
lar effects due to PVC and the polarizations induced by 
deformation and TO mean fields as in the case of ^^Ni. 
Their details will be discussed below. 

Let consider major conclusions emerging from these 
calculations. They are related to the role of polariza- 
tions effects due to deformation and TO mean fields and 
the impact of PVC on the accuracy of the description of 
single-particle spectra, shell gaps, pseudospin doublets 
and spin-orbit splittings. 



Polarization effects due to deformation and 
time-odd mean fields 



FIG. 6: (Color online) Combined polarization effects due 
to deformation and time-odd mean fields. For each single- 
particle state of Figs.[2]l5l the absolute values of the difference 
of the energies obtained in spherical (the Esph energies; column 
"sph" in Figs. [2J5]) and deformed (the Edef+TO energies of the 
calculations with TO mean fields included; column "def-f-TO" 
in Figs.[2]0 calculations are shown. The average polarization 
effects per state Ask {k = tt, i/) are also shown. 



not clear. However, additional bindings due to TO mean 
fields in the nuclei under study obtained in our CRMF 
calculations are within the typical ranges obtained in the 
systematic study of the impact of TO mean fields on the 
binding energies of odd mass nuclei (Ref. Q). In ad- 
dition, the CRMF code employed here has been exten- 
sively and successfully used in the description of rotating 
systems in which TO fields have very large impact on 
the moments of inertia (see Refs. jH, |4l| and references 
quoted therein). These two facts strongly suggest that 
the impact of TO mean fields is correctly described in 
the current work. 

The dominant levels, i. e. the levels with the largest 
spectroscopic factors, as obtained in spherical PVC cal- 
culations, are shown in column "sph-j-PVC" . One can see 
two effects of particle- vibration coupling: (i) the general 
compression of the spectra leading to a better agreement 
with experiment and (ii) in some cases the change of 
the level sequences. With few exceptions, the spectro- 
scopic factors of the dominant levels in the vicinity of 
the Fermi level vary between 0.5 and 0.9 (see Tables Hill 
and IIVI below), thus these states retain basically their 
single-particle nature. 

The impact of particle-vibration coupling on specific 



Combined polarization effects due to deformation and 
TO mean fields are shown in Fig. IHl Two features are 
clearly seen. 

First, combined polarization effects due to deforma- 
tion and time-odd mean fields decrease with an increase 
of mass number. Indeed, with an increase of the size 
of nucleus the relative role of each single-particle orbital 
becomes smaller which leads to the fact that the cal- 
culated quadrupole deformations of odd-mass nuclei de- 
crease with an increase of mass. For example, the aver- 
age proton and neutron deformations of calculated single- 
particle states in ^^Cu and '^''Ni are 0.045 and 0.060, re- 
spectively. On the contrary, these deformations are sig- 
nificantly smaller in ^'^^Bi and ^'^^Sn; they are only 0.0055 
and 0.006, respectively. Thus, odd nuclei neighboring to 
208pj^ remain basically spherical. The decrease of cal- 
culated deformations of odd mass nuclei neighboring to 
doubly magic nuclei with increasing mass clearly indi- 
cates that deformation polarization effects also decrease 
with mass number. In addition, additional binding due 
to TO mean fields decreases with mass (see Ref. [g for 
details). The combination of these two effects is respon- 
sible for the features seen in Fig. IHl 

Fig. ini also reveals that combined polarization effects 
due to deformation and TO mean fields are smaller for 
protons than for neutrons by a factor of approximately 
two. This is especially true for the N = Z ^^Ni and ^™Sn 
nuclei which have the same proton and neutron single- 
particle states. The part of this difference comes from 
the fact that additional binding due to TO mean fields 
is smaller for odd-proton nuclei as compared with odd- 
neutron ones (see Ref. Q). The analysis of Ref. [5] clearly 
indicates that the contributions of the Coulomb force to 
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FIG. 7: (Color online) The deviations of calculated energies 
of the single-particle states from experimental ones obtained 
in the "def+TO" and "hybrid" calculational schemes for in- 
dicated nuclei. The results for proton and neutron states are 
given by solid and open circles. 



the proton single-particle energies in the presence of TO 
mean fields is responsible for this feature. Deformation 
polarization effects show the same features, namely, they 
are smaller for odd-proton nuclei as compared with odd- 
neutron ones. For example, average deformation polar- 
izations per single-particle state are 0.157 MeV and 0.326 
MeV for proton and neutron single-particle states of ^^Ni 
and 0.060 MeV and 0.164 MeV for proton and neutron 
single-particle states of ^°°Sn. Similar to the case of 
polarizations due to TO mean fields, the origin of this 
difference has to be traced back to the Coulomb force 
since compared sets of proton and neutron single-particle 
states are the same in a given N = Z nucleus. 

The majority of the comparisons of the energies of the 
single-particle states obtained in different particle vibra- 
tion coupling models with experiment neglect deforma- 
tion and TO mean field polarization effects induced by 
the odd particle (see, for example, Refs. [H, US Ull ) • The 
current analysis indicates that within the framework of 
CDFT this neglect is more or less justified only for heavy 
nuclei, and it is more justified for proton subsystem than 
for neutron one. 



C. The accuracy of the description of 
single-particle spectra 

The analysis of Figs. [H HI and[5]reveals that the inclu- 
sion of particle vibration coupling substantially improves 
the accuracy of the description of single-particle spectra. 
This statement is quantified in Fig. [7] and Table ID 

Fig. [7] displays the distribution of the deviations be- 
tween calculated and experimental energies of the single- 
particle states in proton and neutron subsystems of ^®Ni, 
^■^•^Sn and ^°^Pb. It clearly shows that in average the 
inclusion of particle- vibration coupling substantially im- 
proves the description of single-particle spectra, and that 
neutron single-particle states are better described than 
proton ones. However, even in the "hybrid" calculations 
there are few (mostly proton) states which deviate from 
experiment by more than 1 MeV. Note that the inclusion 
of particle-vibration coupling can make the agreement 
between theory and experiment worse for some states. 
This is seen, for example, in the case of proton 3si/2 and 
lhg/2 subshells in ^°^Ph (Fig.©. 



TABLE I: Average deviations per state Ae between calculated 
and experimental energies of the single-particle states for a 
proton (neutron) subsystem of a given nucleus. The results 
obtained in the "def-fTO" and "hybrid" calculational schemes 
are shown. 



Nucleus /subsystem 


Aede/+TO [MeV] 


A£hy6„d [MeV] 


"""Ni/proton 


0.76 


0.77 


^®Ni/neutron 


0.89 


0.71 


^^^Sn/proton 


1.02 


0.68 


^''^Sn/neutron 


0.89 


0.39 


^°*Pb/proton 


1.53 


0.84 


^"'^Pb/neutron 


1.00 


0.47 



Average deviations per state Ae between calculated 
and experimental energies of the single-particle states are 
shown in Table H] They are defined as 



_expi 



N 



(5) 



where TV is the number of the states with known experi- 
mental single-particle energies, and sf^ (^i^^) are calcu- 
lated (experimental) energies of the single-particle states. 
One can see that the inclusion of PVC substantially 
improves the description of the single-particle states in 
"'^^^Sn and ^"®Pb. On the contrary, PVC introduces no 
(small) improvement in the description of the proton 
(neutron) single-particle states of ^^Ni. It is interest- 
ing to mention that the best agreement with experiment 
is obtained in the "hybrid" calculations of neutron-rich 
"2Sn. 

While the accuracy of the description of proton and 
neutron states is comparable in ^®Ni, neutron states are 
appreciable better described than proton ones in the "hy- 
brid" calculations of ^■^^Sn and ^^^Pb. The analysis of 
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TABLE II: Extrapolated "experimental" single-particle ener- 
gies of spherical proton and neutron subshells in ^"^"Sn. The 
results of the extrapolations of Refs. [SQ . i65l] are shown. 



Subshell 


e,(MeV) [56] 


£i(MeV) [65] 


Neutrons 


1^*11/2 


-8.6(5) 


-7.8(8) 


2^3/2 


-9.2(5) 


-8.8(8) 


3Sl/2 


-9.3(5) 


-9.3(9) 


2ff7/2 


-10.93(20) 


-10.4(10) 


24/2 


-11.13(20) 


-11.1(10) 


lff9/2 


-17.93(30) 




2pi/2 


-18.38(20) 




Protons 


lff7/2 


+3.90(15) 


2.6(3) 


24/2 


+3.00(80) 


2.8(3) 


lS9/2 


-2.92(20) 


-2.9(3) 


2pi/2 


-3.53(20) 


-3.5(3) 


2p3/2 


-6.38 




1/5/2 


-8.71 





56,,. 100„ 132 _ 208_,, 

Ni Sn Sn Pb 



theses) with the values reported in Ref. [5g, except for 
the energy of the 137/2 proton state. However, estima- 
tion errors are significant which makes detailed compar- 
ison between theory and experiment meaningless. On 
the other hand, one can clearly see in Fig. [3] that PVC 
improves the agreement with experimental estimates. 



FIG. 8: (Color online) The size of proton and neutron shell 
gaps. The results of the "def+TO" and "hybrid" calculations 
are compared with experimental values shown by horizontal 
lines. In the case of ^"''Sn, the experimental shell gaps are 
extracted from extrapolated single-particle energies of Refs. 
[H] (A) and [H (B) [see also Table [TI]. Note that the size 
of the Z = 82 shell gap in ^"^^Pb is almost the same in both 
calculations. 



Figs. S] and [5] strongly suggests that the increases of the 
depth of proton potential in these two nuclei by few hun- 
dreds keV (and as a result the lowering of the energies of 
proton single-particle states by similar energy) will im- 
prove the agreement with experiment. This clearly in- 
dicates that the inclusion of the single-particle energies 
into the fit of the CDFT parametrizations can provide an 
extra information on the depth of the proton and neutron 
potentials and their evolution with particle numbers. 

The detailed comparison between theory and experi- 
ment has not been performed for ^"^Sn since very little 
experimental information is available on neighboring to 
-'^™Sn nuclei with one nucleon less (more) (see, for exam- 
ple, Refs. m, [63). Unfortunately, the extrapolations of 
the single-particle energies from other odd nuclei towards 
^''"Sn do not bring unique values. This is illustrated in 
Tab le llll where the results of the extrapolations of Refs. 
[IE l65| are shown. Note that the estimates of the en- 
ergies of the single-particle states of Ref. coincide 
within the estimation errors (given in Table in paren- 



D. The impact of particle-vibration coupling on 
shell gaps 

The size of the shell gap (between last occupied and 
first unoccupied states in even-even doubly magic nu- 
cleus) is defined as 

<5(Z, N) = miniie,)"'"""'} - maa;{(e, (6) 

where (si)"''"'"^ and {ei)'"^''°^' stand for the energies of the 
single-particle states above and below the shell gap. This 
definition of the shell gap incorporates the information 
from three different nuclei because the single-particle en- 
ergies are defined according to Eqs. ^ and ([2]). The 
advantage of this definition is that it allows to incorpo- 
rate the polarization (due to deformation and TO mean 
fields) and PVC effects on the size of shell gap. 

Fig. [8] compares experimental shell gaps with the shell 
gaps obtained in the "def+TO" and "hy brid" calcula- 
tions. The extrapolations of Refs. [s^, ll^ for the single- 
particle energies in ^'^'^Sn lead to quite different values 
of the proton Z ~ 50 shell gap. It is reasonable to ex- 
pect the same situation for the neutron = 50 shell 
gap, the size of which, however, cannot be determined in 
the extrapolations of Ref. [1^ since no estimates for the 
single-particle states below it are provided (see Table 

The shell gaps are largest in the spherical mean field 
calculations. The gaps become smaller with the inclu- 
sion of each additional type of correlations. The inclusion 
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TABLE III: Spectroscopic factors S of the dominant single-particle levels in odd nuclei surrounding Pb and Sn calculated 
within relativistic particle-vibration coupling model compared to experimental data. The experimental data are taken from 
Refs. [ell (f^Bi, (^He,d) reaction), [SQ] (^^^Bi, the (a,t) reaction), [H (^°'^TL the (d,^He) reaction), M C°''T\, the (t,a) 
reaction), [73] (^"^Pb, the (d,t) reaction), [tJ (^°''Ph, the (^He,a) reaction), ^ (^°^Ph, the (d,p) reaction), [zl (^o^pb, the 
(a, ^He) reaction) and [6l[ (^^^Sn, the (d,p) reaction). For comparison the experimental data from two different reactions are 
presented for odd mass nuclei neighboring to '^"^Pb. 



Nucleus 



'Pb 



'Bi 



'Pb 



'Tl 



State 

2g9/2 

liii/2 

ljl5/2 
3d5/2 
4Si/2 
2g7/2 
3d3/2 



Ih, 



2f 



9/2 



li 



7/2 



13/2 



2f5/2 
3p3/2 
3pi/2 



3p 



1/2 



2f5/2 
3p3/2 



li 



2f 



13/2 



7/2 



lh9/2 

3Sl/2 
2d3/2 

lhii/2 

2d5/2 
lg7/2 



t>th 

0.85 
0.89 
0.66 
0.89 
0.92 
0.87 
0.89 



0.90 
0.87 
0.86 
0.82 
0.64 
0.38 

0.84 
0.86 
0.80 
0.68 
0.22 



0.78±0.1 
0.96±0.2 
0.53±0.2 
0.88±0.1 
0.88±0.1 
0.78±0.1 
0.88±0.1 



0.88 


1.17 


[68 


0.80 


0.78 


0.78 


[68 


0.76 


0.63 


0.56 


[68 


0.74 


0.61 


0.88 


[68 


0.57 


0.62 


0.67 


[68 


0.44 


0.37 


0.49 


[68 


0.20 



1.13 
1.00 
1.04 
0.89 



0.95 
1.15 
0.89 
0.62 
0.40 



M 

M 

69] 

69] 
69] 



^exp 

0.94 [73] 
1.05 [73] 
0.57 [73] 



1.08 [71 
1.05 [71 
0.95 [71 



0.61 
0.64 



0.85 [35 
0.90 [35 
0.88 [3^ 
0.63 [3^ 
0.27 [35 



Nucleus 



'Sn 



'Sb 



^Sn 



^In 



State 



2f7/2 
3p3/2 
lh9/2 

3pi/2 

2f5/2 



2d3/2 
lllll/2 

3Si/2 



2di 
h 



5/2 
i7/2 



ig9/2 
2p3/2 

2pi/2 

lf5/2 



0.89 
0.91 
0.88 
0.91 
0.89 



lg7/2 


0.86 


2d5/2 


0.82 


2d3/2 


0.63 


lhii/2 


0.79 


3Si/2 


0.70 



0.88 
0.86 
0.87 
0.70 
0.72 



0.85 
0.70 
0.85 
0.37 



^exp 

0.86±0.16 
0.92±0.18 

1.1±0.3 
1.1±0.2 



of polarization effects due to deformation and TO mean 
fields decreases slightly tfieir sizes (see Figs. [U [31 S] and 
[5]). Thus, we compare with experiment in Fig. [8] only the 
PVC ( "hybrid" ) and best mean field ( "def-|-TO" ) calcula- 
tions. This figure shows that particle-vibration coupling 
decreases substantially the size of the shell gaps. The 
effect is most pronounced in ^^Ni where PVC decreases 
the proton Z = 28 and neutron = 28 shell gaps by 
almost 2 MeV. On the contrary, the effect of particle- 
vibration coupling is least pronounced in ^"^Pb in which 
it decreases the size of the neutron iV = 126 gap by only 
few hundreds keV and has almost no impact on the pro- 
ton Z = 82 shell gap. 

Comparing experimental and calculated shell gaps in 
^^Ni, ^■^^Sn and ^"*Pb, one can conclude that both calcu- 
lational schemes provide similar accuracy of the descrip- 
tion of experimental data. However, experimental shell 
gaps are typically overestimated in the "def-|-TO" calcu- 
lations. The exceptions are the Z = 50 gap in ^^^Sn and 
Z = 82 gap in ^"^^Pb. On the contrary, with exception 
of the N = 126 gap in ^°^Pb, the size of experimental 
gaps is underestimated in the calculations with particle- 



vibration coupling. 



E. Spectroscopic factors 

The calculated spectroscopic factors for odd nuclei 
with doubly-magic cores are compared with experiment 
in Tables Hill and ITVl In odd nuclei neighboring to ^"^^Pb, 
the experimental spectroscopic factors are reasonably 
well reproduced in the PVC calculations. Both in ex- 
periment and in calculations the proton states are found 
to be somewhat more fragmented than the neutron ones. 
When comparing theory with experiment one should keep 
in mind that the experimental spectroscopic factors de- 
pend considerably on the parameters used in the model 
analysis and on employed reaction, which is clearly seen 
in the staggering of experimental values. Note that the 
results obtained for ^°*Pb with the NL3* parametriza- 
tion are very close to those obtained previously in Ref. 
[Hi with the NL3 force. 

For other nuclei, the experimental data are available 
only for ^^^Sn (Table HIJ) and ^^Ni (Table |IV|. The level 
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TABLE IV: Spectroscopic factors S of the dominant single-particle levels in odd nuclei surrounding ^""Sn and ^^Ni calculated 
within relativistic particle- vibration coupling model. They are compared to the experimental values extracted by means of the 
{d,p) transfer [sj and one-neutron knockout [t^I reactions in the case of ^^Ni. 



Nucleus 


State 


Sth 


Nucleus 


State 


Sth 


Sexp 34J 


Se:rp 174] 


"'^Sn 


2d5/2 
lg7/2 
2d3/2 
3Si/2 

lhii/2 


0.85 
0.85 
0.78 
0.81 
0.80 


5^Ni 


2p3/2 
lf5/2 

2pi/2 

lg9/2 


0.83 
0.79 
0.76 
0.79 


0.95±0.29 
1.40±0.42 
1.00±0.30 


0.58±0.11 




2d5/2 
lg7/2 
2d3/2 
3Si/2 

lhii/2 


0.87 
0.86 
0.83 
0.87 
0.81 


^■^Cu 


2p3/2 
lf5/2 

2pi/2 

lg9/2 


0.85 
0.80 
0.80 
0.80 






9«Sn 


lg9/2 
2pi/2 
2P3/2 
lf5/2 
lf7/2 


0.84 
0.85 
0.71 
0.62 
0.11 


5^Ni 


lf7/2 
2Si/2 
ld3/2 
ld5/2 


0.78 
0.71 
0.62 
0.20 






*In 


lg9/2 

2pi/2 

2p3/2 
11,5/2 
lf7/2 


0.85 
0.86 
0.73 
0.66 
0.14 


Co 


lf7/2 
2Si/2 
ld3/2 
ld5/2 


0.78 
0.73 
0.64 
0.20 
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FIG. 9; (Color online) Spin-orbit splitting energies AEu of 
the spin-orbit doublets in ^"^^Pb. Insert shows the colors used 
for different values of orbital angular momentum I. Theo- 
retical results are shown by symbols, while the experimental 
values by solid lines. Solid (open) symbols are used for spin- 
orbit doublets which are build from the particle (hole) states 
with respect of the ^"^Pb core. The stars and crosses are 
used for the spin-orbit doublets which involve both particle 
and hole states. 



of agreement between theory and experiment in Sn is 
comparable with the one seen in odd nuclei neighbouring 
to 208pb. The situation in "Ni is more controversial. For 
example, the sum rule ^ is strongly violated in the case 



of the 1/5/2 state in experimental data of Ref. [t^. In 
addition, experimental spectroscopic factors for the 2p^/2 
state in ^^Ni extracted by means of the {d,p) transfer [SJ] 
and one-neutron knockout [t^ reactions differ by a factor 
of almost two. Currently there is no satisfactory expla- 
nation for this difference (Ref. [l^]). The results of the 
PVC calculations are lower than the experimental values 
obtained in the (d, p) transfer reaction but higher than 
the ones obtained in the one-neutron knockout reaction. 
On the other hand, the spectroscopic factors obtained in 
the PVC calculations are very close to the ones obtained 
in the spherical shell model calculations in full fp model 
space which employ GXPFIA interaction (see Table II 
in Ref. [34]). Note that the latter calculations provide 
best agreement with the data in the Ni isotopes [3^. 
One should also mention that non-relativistic particle- 
vibration calculations of Ref. [zl] give the spectroscopic 
factors for the 2^3/2, I/5/21 and 2pi/2 states which are 
very close to 0.6. 



F. The impact of particle-vibrational coupling on 
the spin-orbit splittings 

Covariant density functional theory naturally describes 
the spin-orbit interaction in nuclei [l|, 01 , which is a rela- 
tivistic effect that has to be added phenomenologically in 
non-relativistic models. The comparison of experimental 
spin-orbit splittings with calculations has been performed 
only on the mean field level so far (see, for example, Refs. 
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FIG. 10: (Color online) The same as in Fig. [9] but for the 
spin-orbit doublets of ^^^Sn. 



[1,[T^). However, the current work clearly indicates that 
particle-vibration coupling has an impact on the ener- 
gies of the single-particle states. Thus, it is reasonable to 
expect that it will also modify the spin-orbit splittings. 

In order to understand the impact of polarization 
effects (due to deformation and TO mean fields) and 
particle-vibration coupling on spin-orbit splittings, the 
results of the "sph", "sph-fPVC", "def", "def-fTO", and 
"hybrid" calculations are compared with available exper- 
imental data on the spin-orbit doublets in Figs. [^[TUl and 

HH 

The impact of deformation and TO mean field polar- 
ization effects on spin-orbit splittings is rather modest in 
208pj-, (Pig, since these effects are smallest in heavy 
systems (see Sec. IIVB[) . In addition, both of them lead 
to a compression of the single-particle spectra. As a con- 
sequence, they act in the same direction for the members 
of the spin-orbit doublet which is built either from hole 
(below the shell gap) or particle (above the shell gap) 
states. Moreover, these polarization effects are more or 
less similar in magnitude for the members of such spin- 
orbit doublets. As a consequence, they cancel to a large 
extent when spin-orbit splitting is calculated as the differ- 
ence of the energies of the members of the spin-orbit dou- 
blet. Such behavior is typical for the spin-orbit doublets 
based on the orbitals with low orbital angular momen- 
tum both members of which arc located either below 
or above the shell gap. For this type of the spin-orbit 
doublets, the impact of polarization effects on spin-orbit 
splittings is also small in -'^'^^Sn (Fig. [TU)) . However, it 
becomes more pronounced (especially, in neutron sub- 
system) in s^Ni (Fig.lTTt. 

The situation is different for the spin-orbit doublets 
built on the orbitals with high orbital angular momentum 
I. For such doublets, one member is located above the 
shell gap, while another below the shell gap. The com- 
pression of the single-particle spectra due to deformation 
and TO mean field polarization effects leads to the de- 



crease of spin-orbit splitting. This decrease is especially 
pronounced for light nuclei and for neutron subsystem in 
a given nucleus (compare Figs. [TTl [TUl andlH]) because of 
particle number and subsystem (proton or neutron) de- 
pendencies of polarization effects discussed in Sec. IIVBI 
The impact of particle- vibration coupling on spin-orbit 
splittings is state-dependent. It is small or modest in the 
proton / — 3 (2/7/2 — 2/5/2) and neutron / = 1 — 
3pi/2), I = 2 (2^5/2 - 2^3/2), and / = ^ (2^9/2 - 2^7/2) 
doublets of ^o^pb (Fig. proton / = 2 (2^5/3 - 2^3/2) 
and neutron / = 1 (87)3/2 - 'ipi/2), I = 3 (2/7/2 - 2/5/2) 
doublets of ^"^^Sn (Fig. [TO)) and neutron and proton / = 
1 (2p3/2 - 2pi/2) doublets in ^^Ni (Fig. HI]). In these 
doublets, the corrections to the energies of the members 
of the spin-orbit doublet due to PVC are more or less the 
same. 

However, this is not always the case for the doublets 
based on the orbitals with low orbital angular momentum 
I. For example, the corrections due to PVC differ signif- 
icantly for the proton 3^3/2 and 3c?5/2 states of ^°^Pb 



(Fig. [S]). As a consequence, particle- vibration coupling 
has significant effect on the spin-orbit splitting of the 
proton I — 2 (3c?5/2 — 8^3/2) doublet. In addition, spin- 
orbit splittings of the proton / = 1 (8^3/2 — 8^1/2) and 
neutron I = 3 (2/7/2 - 2/5/2) doublets in ^ospb (Fig. 
[9]) as well as of neutron / — 2 (2c?5/2 — 2^3/2) doublet in 
^■^^Sn (Fig.fTU)) are strongly affected by particle-vibration 
coupling. 

The effect of particle-vibration coupling is especially 
pronounced for the spin-orbit doublets with high orbital 
angular momentum I, the one member of which is lo- 
cated below the shell gap and another above the shell 
gap. These are the proton I ~ 5 (l/iii/2 — 1^9/2) ^md 
neutron I = 6 (IH3/2 — 1*11/2) doublets in ^°^Pb (Fig. 
[H]), the proton I — 4 (I59/2 ~ I57/2) ^nd neutron / = 5 
(l/iii/2 - lft-9/2) doublets in i32Sn (Fig. nop and the pro- 
ton and neutron I = 3 (I/7/2 — I/5/2) doublets in ^^Ni 
(Fig. [TT]). For these doublets, the spin-orbit splittings 
are largest in spherical calculations. The inclusion of po- 
larization effects due to deformation and TO mean fields 
decreases these splittings. The PVC coupling further re- 
duces the energy splittings by approximately 1 MeV in 
medium and heavy mass nuclei (Figs. [9l and [TO|) and by 
approximately 1.5 MeV in ^^Ni fFig. [TT|) . 

Fig. [T2] presents percentual deviations S{AEis) of the 
calculated spin-orbit splittings from experimental ones 
defined as 



ae: 



■exp 



X 100% 



(7) 



where AEis is the spin-orbit splitting energy. Negative 
(positive) values of 5{AEis) indicate that experimental 
spin-orbit splitting is underestimated (overestimated) in 
the calculations. One can see that in average the mean- 
field results are more or less evenly scattered around zero 
percentual deviation. On the contrary, with few excep- 
tions particle-vibration coupling calculations underesti- 
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mate experimental spin-orbit splittings. Spin-orbit dou- 
blets built on the orbitals with high orbital momentum 
I, involving the members across the shell gap, are rather 
well (within the 20% deviation) described in both calcu- 
lations. These doublets are characterized by large spin- 
orbit splittings AEis in the range of 5-7 MeV, so the 20% 
deviation means that the absolute deviation from experi- 
ment is typically less than 1 MeV. The largest percentual 
deviations are observed for the spin-orbit doublets built 
on the orbitals with low orbital angular momentum /. 
These doublets are characterized by low spin-orbit split- 
tings so relatively low absolute deviations from experi- 
mental values of the order of few hundreds keV result in 
appreciable percentual deviations. 

The analysis of the results of calculations allows to 
conclude that the inclusion of particle-vibration coupling 
decreases the accuracy of the description of spin-orbit 
splittings. This is clearly visible when absolute and/or 
percentual deviations per doublet are compared in the 
mean field and PVC calculations. The absolute devia- 
tions per doublet are 0.34 MeV [0.50 MeV], 0.23 MeV 
[0.56 MeV] and 0.26 MeV [0.45 MeV] in the mean field 
("def-|-TO") [particle- vibration coupling ("hybrid")] cal- 
culations in ^^Ni, ^'^^Sn and ^"^Pb, respectively. The per- 
centual deviations per doublet are 11.8 % [10.3%], 14% 
[21.5%] and 19.3% [36.4%] in the mean field ("def-fTO") 
[particle- vibration coupling ("hybrid")] calculations in 
^^^Ni, i^^Sn and ^ospb, respectively. 

These results are not surprising and should not be 
viewed negatively. First of all, as discussed in Sees. 
and IIV El the experimental levels are not pure single- 
particle levels, and as such they are closer in nature to 
the calculated levels of the PVC model. Second, rela- 
tivistic description 01 predicts a definite connection be- 
tween the Dirac effective mass m^/m of the nucleon in 
the kinetic energy and the strength of the spin-orbit force 
because the same effective mass appears in both terms. 
In addition, the spin-orbit term is sensitive to the spa- 
tial variations of the effective mass. Particle- vibration 
coupling affects the effective mass of the nucleon and as 
a result it has an impact on calculated spin-orbit split- 
tings. Third, these results point to a new direction of 
improving the covariant energy density functionals. This 
is because empirical energy spacings between spin-orbit 
partner states in finite nuclei determine a relatively nar- 
row interval of allowed values for the Dirac effective mass 
0.57 <m}j/m< 0.61 on the mean field level [Hi. The 
observed impact of particle- vibration coupling on spin- 
orbit splittings suggests that this interval of the Dirac 
effective mass m^/m may in reality be broader. 



G. The impact of particle-vibtational coupling on 
the pseudospin doublets 
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FIG. 11: (Color online) The same as in Fig. [9] but for the 
spin-orbit doublets of ^®Ni. 
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FIG. 12: (Color online) Percentual deviations of the spin- 
orbit splittings in the "def-|-TO" (red open circles connected 
by dashed lines) and "hybrid" (black solid circles connected 
by solid lines) from the experimental values. 



of attractive scalar S and repulsive vector V relativis- 
tic mean fields in which the nucleons move'^. Pseu- 
dospin doublets have nonrelativistic quantum numbers 
{nr,l,j = 1 + 1/2) and {ur - 1,1 + 2, j = I + 1/2) where 



It was shown in Ref. 6] that quasidegenerate pseu- 
dospin doublets in nuclei, discovered more than 40 years 
ago (76l . [ttI , arise from the near equality in magnitude 



Extensive review of the manifestations of pseudospin symmetry 
in different physical systems is presented in Ref. [7S| . 
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rir, I and j are the single-nucleon radial, orbital and total 
angular momentum quantum numbers, respectively^. 

The pseudospin doublet can be characterized by pseu- 
dospin doublet splitting energy which is defined as the 
difference of the energies of pseudospin doublet members 

= E„^jj=i+i/2 — E(^„,^_l),{l+2),j=l+3/2- (8) 

With this definition of AE, the low-/ (high-/) member of 
the pseudospin doublet is higher in energy than its high- 
/ (low-/) counterpart ii AE > {AE < 0). Note that 
the calculations fail to reproduce the relative order (in 
energy) of the members of pseudospin doublet if the signs 
of AE are different in the calculations and experiment. 

Fig. [13] compares available experimental data on pseu- 
dospin doublet splitting energies in ^^Ni, ^^^^Sn and ^°^Pb 
with PVC ("hybrid") and best mean field ("def-^TO") 
calculations. Note that the AE values obtained in the 
"sph" and "def" calculations deviate from the ones ob- 
tained in the "def+TO" calculations by less than 200 keV 



* Note that in the current manuscript the single-particle subshells 
are labelled by the labels which use (rir -I- 1) in the first position 
of the label. 



(in the majority of the cases the difference is less than 
100 keV). This is because both members of pseudospin 
doublet are located either above or below the shell gaps. 
The polarization effects due to deformation and TO mean 
fields act in the same direction for these states. As a 
result, these effects cancel each other to a large degree 
when the difference of the energies of the members of 
pseudospin doublet is taken. 

Fig. [13] clearly shows that particle- vibration coupling 
substantially improves the description of splitting ener- 
gies in pseudospin doublets; the average deviations from 
experiment for AE are 0.93 MeV and 0.6 MeV in the 
mean field and PVC calculations, respectively. There 
are still some doublets in which the energy splitting AE 
is poorly reproduced in model calculations. These are 
7r2d5/2lg7/2 and 7r2/7/2l/i9/2 doublets in ^°^Pb (see Figs. 
[5] and [T3]) , the AE values of which deviate from experi- 
ment by 2.18 MeV (1.66 MeV) and 2.16 MeV (1.2 MeV) 
in the mean field (PVC) calculations, respectively. We 
did not find clear explanation for such large differences. 

Proton and neutron pseudospin doublets with the same 
single-particle structure have similar splitting energies 
in experiment (see pseudospin doublets indicated by the 
brackets with arrows in Fig. 1131) . This feature is rather 
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well reproduced both in the mean field and PVC cal- 
culations for the proton and neutron 2P3/2I/5/2 pseu- 
dospin doublets of ^^Ni. On the contrary, the mean 
field ("def+TO") calculations completely fail to repro- 
duce this feature for the proton and neutron 2^5/21.97/2 
pseudospin doublets in ^■^^Sn and 2/7/21/19/2 pseudospin 
doublets in 208p|-, (Fig.[T31); the calculated splitting ener- 
gies for neutron doublets are by more than 1 MeV larger 
than the ones for proton doublets. However, in the PVC 
calculations the splitting energies for these pairs of pro- 
ton and neutron pseudospin doublets are similar in agree- 
ment with experiment. This is a consequence of the fact 
that particle-vibration coupling decreases (as compared 
with mean field calculations) the splitting energies in neu- 
tron pseudospin doublets by more than 1 MeV leaving at 
the same time the splitting energies in proton pseudospin 
doublets unchanged. 



H. The 120 nucleus - the center of the island of 
stabihty of superheavy nuclei 

The superheavy ^^^120 nucleus has been predicted to 
have proton Z = 120 and neutron = 172 spherical shell 
closures within the RMF theory Thus, it represents 
the center of the "island of stability" of shell-stabilized 
superheavy nuclei. The analysis of the deformed one- 
quasiparticle states in the A ^ 250 mass regions in Ref. 
[23 | supports the presence of the Z = 120 shell gap in 
superheavy nuclei. In addition, it indicates N = 172 as a 
likely candidate for magic neutron number in superheavy 
nuclei. 

This nucleus represents a challenge for future experi- 
mental synthesis since it is located at the limits of ac- 
cessibility by available cold fusion reactions. Therefore, 
as accurate as possible estimations of its characteristics 
are needed from the theoretical side. First results for 
its single-particle spectra obtained within the relativistic 
particle-vibration coupling model are presented in Fig. 
1141 They are compared with the results obtained on 
the mean field level. Note that we restrict ourselves in 
the case of the ^^^120 nucleus to spherical calculations. 
This is because the polarization efi^ects due to deforma- 
tion and TO mean fields decrease with mass (see Sect. 
IIVB|) . and thus their impact on the single-particle spec- 
tra is expected to be rather small in the ^^^120 nucleus. 

The density of the single-particle states is substantially 
larger than in lighter nuclei at the mean field level. As a 
consequence, this nucleus exhibits a very rich spectrum 
of low-lying collective phonons already within the RRPA. 
For example, the lowest vibrational 2f mode appears at 
1.41 MeV with B(E2)t = 7.1 xlO^ e^ fm^, accompa- 
nied by several rather collective modes at 3.18 MeV, 5.25 
MeV, 6.74 MeV etc. Similar picture is obtained for other 
multipolarities. Although contributions to the nucleonic 
self-energy from the high-lying modes decrease quickly 
with energy, we have included phonons with all J'^ val- 
ues mentioned in Sec. IIII Dl below 15 MeV, altogether 



about 100 phonons. 

Fig.ll4lshows that the particle- vibration coupling gives 
rise to the general compression of the spectra. As a con- 
sequence, the size of the Z ~ 120 shell gap decreases 
from 3.35 MeV (in mean field calculations) down to 2.67 
MeV (in the PVC calculations), and the size of the neu- 
tron shell gap at TV = 172 from 2.42 MeV down to 1.83 
MeV. Although these gaps are smaller by 20-25% than 
the ones obtained in the mean field calculations, they are 
still significant. As a result, the ^^^120 nucleus still re- 
mains doubly magic nucleus representing the center of 
the "island of stability" of superheavy nuclei even when 
the correlations beyond mean field are taken into account. 

V. CONCLUSIONS 

The relativistic particle-vibration model has system- 
atically been applied in combination with cranked rela- 
tivistic mean field approach for the study of the impact 
of surface vibrations on single-particle motion. The po- 
larization effects in odd mass nuclei due to deformation 
and time-odd mean fields have been treated in the CRMF 
framework. The main results can be summarized as fol- 
lows: 

• Particle-vibration coupling has to be taken into ac- 
count when model calculations are compared with 
experiment since this coupling is responsible for ob- 
served fragmentation of experimental levels. The 
inclusion of particle-vibration coupling substan- 
tially improves the description of the energies of 
dominant single-particle states in i32Sn and ^ospb. 
However, the accuracy of the description of single- 
particle spectra in ^^Ni is similar in the mean field 
and PVC calculations. Note that dominant neutron 
single-particle states are in average better described 
than proton ones in the PVC calculations. 

• The polarization effects in odd-mass nuclei due 
to deformation and time-odd mean fields induced 
by odd particle are important. They have to be 
taken into account when experimental and calcu- 
lated single-particle energies are compared. How- 
ever, they are usually neglected when the results of 
particle-vibration coupling model calculations are 
compared with experiment. The current analysis 
indicates that within the framework of CDFT this 
neglect is more or less justified only for heavy nu- 
clei, and it is more justified for proton subsystem 
than for neutron one. 

• Particle-vibration coupling leads to a shrinkage of 
the shell gaps. The size of the shell gaps is typically 
underestimated in the PVC calculations as com- 
pared with experiment and overestimated in the 
mean field calculations. 

• The inclusion of particle-vibration coupling de- 
creases the accuracy of the description of spin-orbit 
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FIG. 14: (Color online) Single-particle spectra of the ^"^120 nucleus. Column 'sph' shows the single-particle spectra obtained in 
spherical RMF calculations. Column 'sph-|-PVC' shows the spectra obtained in spherical calculations within the PVC model. 



splittings. The analysis suggests that the mean 
field constraint on the allowed range of the Dirac ef- 
fective mass 0.57 < mj,/m < 0.61 is too restrictive 
and the range of the /m values can be broader 
in the models which take into account the correla- 
tions beyond the mean field. 

• Particle-vibration coupling substantially improves 
the description of splitting energies in pseudospin 
doublets as compared with mean field calculations. 
Observed similarity of the splitting energies of pro- 
ton and neutron pseudospin doublets with the same 
single-particle structure in medium and heavy mass 
nuclei can only be reproduced when the particle- 
vibration coupling is taken into account. 

• The spherical shell closures in superheavy nuclei are 
still found at proton Z ~ 120 and neutron N — 172 
numbers even when particle-vibration coupling is 
taken into account. However, the size of these gaps 
becomes smaller (as compared with mean field val- 
ues) in the presence of particle-vibration coupling 
due to general compression of the spectra caused 
by the increase of the effective mass of the nucleon 
at the Fermi level. 



• The NL3* parametrization employed in the cur- 
rent work has been adjusted at the mean-field level. 
The remaining discrepancies between the results 
of particle vibration coupling model and experi- 
ment clearly suggest that this parametrization is 
not completely adequate for the description of the 
energies of the single-particle states. We believe 
that this statement is not limited to NL3''' but it is 
also valid for all existing CDFT parametrizations 
which have been adjusted at the mean-field level. 
This calls for the parametrizations specifically tai- 
lored to describe single-particle degrees of freedom 
in the models taking into account the correlations 
beyond the mean field. 
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